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Abstract 

We describe the implementation of Quantum Electrodynamic (QED) evolution at Leading Order 
(LO) along with Quantum Chromodynamic (QCD) evolution at Next-to-Leading Order (NLO) 
in the CTEQ-TEA Global analysis package. The inelastic contribution to the photon Parton 
Distribution Function (PDF) is described by a two-parameter ansatz, coming from radiation off 
the valence quarks, and based on the CT14 NLO PDFs. Setting the two parameters to be equal 
allows us to completely specify the inelastic photon PDF in terms of the inelastic momentum 
fraction carried by the photon, Pq, at the initial scale Qo = 1-295 GeV. We obtain constraints on 
the photon PDF by comparing with ZEUS data [l0| on the production of isolated photons in deep 
inelastic scattering, ep —)• ey + X. For this comparison we present a new perturbative calculation of 
the process that consistently combines the photon-initiated contribution with the quark-initiated 
contribution. Gomparison with the data allows us to put a constraint at the 90% confidence level 
of Pq < 0.14% for the inelastic photon PDF at the initial scale of Qo = 1-295 GeV in the one- 
parameter radiative ansatz. The resulting inelastic GT14QED PDFs will be made available to the 
public. In addition, we also provide GT14QEDinc PDFs, in which the inclusive photon PDF at 
the scale Qo is defined by the sum of the inelastic photon PDF and the elastic photon distribution 
obtained from the Equivalent Photon Approximation. 
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I. INTRODUCTION 


The high precision of current collider data requires comparable precision in the phe¬ 
nomenological predictions. The state of the art in high-energy calculations is at next-to- 
next-to-leading order (NNLO) in Quantum ChromoDynamics (QCD). Consequently, major 
efforts have been undertaken to produce NNLO parton distribution functions (TDFs) from 
a global analysis of the available data. These include the CT14NNLO PDFs l| as well as 


others all of which include LHC data in the determination of the PDFs. 

In this paper we describe the introduction of QED evolution at leading order (LO) with 
the next-to-leading order (NLO) QCD evolution in the same CTEQ global analysis package 
that was used to produce the CT14 PDFs [l|. Past studies of QED effects in global analysis 


have been done by the MRST and the NNPDF p groups. We have checked our code 


against other QED-I-QCD evolution codes 


and hnd good agreement. 


The MRST and NNPDF analyses used different approaches for modeling the photon 
PDF. The MRST group used a parametrization for the photon PDF based on radiation off 
of “primordial” up and down quarks, with the photon radiation cut off at low scales by 


constituent or current quark masses The NNPDF group used a more general photon 
parametrization, which was then constrained by high-energy W, Z and Drell-Yan data at the 
LHC j^. They found constraints on the size of the photon PDF, which was still consistent 
with zero at the initial scale of \p2 GeV. As discussed by Martin and Ryskin j^, the photon 
PDF has a large elastic contribution in which the proton remains intact, in addition to the 
inelastic contribution in which the proton breaks into a multihadron hnal state.^ Neither 
MRST nor NNPDF addresses these separate contributions to the photon PDF, although 
we can assume that the NNPDF photon is inclusive, containing both inelastic and elastic 
components, since it was constrained using inclusive Drell-Yan and vector boson data. 

Given the limited amount of data to constrain the shape of the photon PDF, we will use 
a generalization of the MRST approach. We parametrize the inelastic contribution to the 
photon at the initial scale^ Qq = 1.295 GeV by 

,n/\\ 


fy/p(^x^ Qo) 2 ^ ^ ^ (^) ^ ^ (^) 


where P-yq o /°(x) is the convolution of the quark-to-photon splitting function P-yq^x) with 
the “primordial” quark distribution f^{x), which we take to be the initial GT14 NLO up 
and down valence distributions. We then set Au = Ad to obtain a single parameter family of 
photon distributions, which we can label by their initial inelastic momentum fraction For 
comparison, in analogy with the MRST approach, we will also show results for a “Gurrent 
Mass” (GM) photon distribution, given by dehning A, = and setting the Qi to 

the quark current masses; i.e., Qu = nxu = 6 MeV and Qd = nvd = 10 MeV. 


^ In Ref. [Ol these two contributions are referred to as “coherent” and “incoherent”, respectively. 

^ The initial scale Qo = 1-295 GeV is the same as that used in the standard CT14 PDF sets, and was chosen 
to be just below the input charm pole mass of rric = 1.3 GeV. 
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We will constrain the inelastic photon PDF using data on Deep Inelastic Scattering (DIS) 


with isolated photons from the ZEUS collaboration 10|. The advantage of using this process 


is that the initial-state photon contributions are at leading order in the perturbation expan¬ 
sion. In contrast, the initial-state photon contribution to Drell-Yan or W and Z production 
is suppressed by factors of (a/ag) relative to the leading quark-antiquark production. How¬ 
ever, to use the DIS-plus-photon data we will Erst need to address some technical issues 
relating to the combination of different subprocess contributions to the observed hnal state. 

The organization of this paper is as follows: In Sec. HIl we describe the inclusion of QED 
evolution in the CTEQ global analysis code and give more details about our initial PDF 
parametrizations. In Sec. IIIII we discuss constraints on the photon PDF coming from the 
CT14 global analysis data set and from the ZEUS DIS with isolated photon data. In this 
section we present a new calculation for the DIS-plus-isolated-photon process, which con¬ 
sistently combines the photon-initiated contribution with the quark-initiated contribution. 
We show that this data gives signihcant constraints on the initial photon PDF. In Sec. HV] 
we discuss our Endings and give conclusions. We also include an Appendix where we show 
comparisons between our QCD-I-QED evolution code and other publicly-available codes. 


II. INCORPORATION OF QED EFFECTS IN CTEQ-TEA GLOBAL ANALYSIS 

In this section we discuss the implementation of the QED evolution and the initial photon 
PDF in the context of the CTEQ-TEA global analysis program. 
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A. QCD-plus-QED evolution 


The evolution of the PDFs, including QED contributions at leading order (LO) 

and QCD contributions at higher orders, is described by the equations: 


dfq. 

dt 


dt 


dh 

dt 

df^ 

dt 




{d^qiQj ° fqj + d^qiQj ° fqj) + d^qg ° fg 


o f + o f 


271 - ° fgj d- Pqiqj ° fqj) + Pqg ° fg 


-4- — e^ I P^°^ O f- + P® O f 


a 


Pgg ° fg P ^sq ° (/■ji + 


( 2 ) 


where t = In fff, the indices i and j run over active quark flavors, and the convolution is 
dehned by 


dz dy6{zy-x)Pab{z)fb{y,IJ^F)- (3) 

Jo 

The QCD splitting functions, given by 


{PabO fb)ix,yF) = 



n 


( 4 ) 
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are known up to n = 2, next-to-next-to-leading order (NNLO) 
splitting functions can be extracted from the LO QCD splitting functions, giving 


ij. The LO QED 



\ z^ 3 . . 

A, + (1 - 

1 +( 1-^)2 


z) , 


( 5 ) 

( 6 ) 


where = 3 is the number of colors. 
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We have modified the Fortran NLO evolution code evolve, which was used for previous 

and CT09 17|), to include the LO QED contri- 


CTEQ-TEA PDFs (CTEQ6-6.6 [I3j-[16 


butions. This code solves the evolution equations directly in x-space, so that the only new 
technical issue introduced by the QED corrections is the separation of the quark singlet 
distributions into separate up and down contributions, based on the quark charges. We have 
checked our evolution code against the public QCD+QED codes, partonevolution BQ 
and APFEL we find good agreement. Details of this comparison are given in the 

Appendix. 


B. Initial Photon PDFs 


The initial photon PDF at the scale Qq is a nonperturbative input that must be obtained 
by a fit to data. Even a choice of zero initial photon PDF is ambiguous, since it depends on 
the arbitrary scale Qq. So far there have been two different approaches to the initial photon 
PDF. In the MRST analysis B fhe initial photon PDF was given by an ansatz, obtained 
from radiation off primordial valence up and down quark distributions, cut off at low scales 
given by the current quark masses, = 6 MeV and rrid = 10 MeV, or by constituent quark 
masses mu = mu = 300 MeV. Alternatively, the NNPDF approach B was to use a general 
parametrization for the initial photon PDF to be constrained by high-energy W, Z, and 
Drell-Yan production at the LHC. 

In this work we will use a generalization of the MRST ansatz, but we must first address a 
subtlety of the photon PDF. Unlike the case for colored partons, the photon PDF has a large 
elastic component, in which the proton remains intact [91]. This is in addition to the inelastic 
component, in which the proton dissociates. The elastic component can be parametrized by 
the Equivalent Photon Approximation (EPA) 19(|, which involves an integration over the 
proton electromagnetic form factors. For this work we focus on the inelastic component, 
which we parametrize by a radiative ansatz, but with free parameters to be fit by data. 
Given the weak constraints from data on the photon PDF, we find it useful to limit the 
number of parameters to one or two for the time being. We shall see that the ZEUS DIS- 
plus-isolated-photon data 1^ constrains the inelastic photon PDF roughly in the range 
10-3 <x<2- 10-2 for 16 < Q2 < 300 GeV^. 

We parametrize the inelastic contribution to the initial photon PDFs in the proton and 
neutron by 


Qo) 2^ ° ^ (^) A ° ^ 

f'r/n{x,Qo) = ■^(^AuelP^god°{x) + AdelP^qOu%x)^ , (7) 


where and are “primordial” valence-type distributions in the proton, and the initial 
photon PDF in the neutron is obtained by an approximate isospin symmetry. Defining 
Ai = we can trade the parameters Aj for mass scales Q*, and we see that the 
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nonperturbative inputs f-y/{p,n){.x,Qo) are modeled by the radiation of a single photon off 
the “primordial” quarks, with a collinear cutoff given by the scales Qi. The MRST initial 
photon PDFs can be obtained from this parametrization by setting Qq = 1 GeV, using the 
functions for uq and do given in Ref. [^, and letting Qu and Qd be either the constituent or 
current quark masses. For our analysis, we use Qo = 1.295 GeV as in CT14, and we set 

uoix) = u^{x,Qo) = fu/pix,Ql) - fu/pix,Ql) 

do{x) = dFy{x,Qo) = fd/p{x,Ql) - fd/p{x,Ql) , ( 8 ) 

the initial up and down CT14 NLO valence distributions in the proton. 

The presence of the photon PDF violates isospin between the neutron and proton. Con¬ 
tinuing with the radiative ansatz and working to hrst order in a, we can neglect the isospin 
violation in the gluon and sea-quark PDFs [51] and use 

Qo) fg/n{,X, Qo) 

fq/pix,Qo) = fq/nix,Qo) for q = u,d,s,c,b (9) 

fq/pix,Qo) = fq/n{x,Qo) for q = s,c,b. 

For the valence quarks at hrst order in a, the radiative ansatz plus approximate isospin 
symmetry implies 

cv ~ 

vFy{x, Qo) ~ vP{x) + —AuelPqg o u°{x) 

Ztt 

dy{x, Qo) ^ d°(x) -F ^AdelPqq O d^{x) 

Ztt 

O/ ~ 

UyQc, Qo) ~ + Tr^uOiPqq o d^{x) (10) 

ZTT 

dy{x,Qo) ~ M°(x) ^Ade\PqqOU^{x) . 

This suggests a consistent set of PDFs for the valence quarks in the neutron given by 

Uy{x, Qo) = dy(x, Qo) + ^ {Auel - AdC^) Pqq o d^{x) 

dv{x, Qo) = My(x, Qo) + ^ {Adel - A^el) Pqq o u^{x) . (11) 

Note that Eqs. o and flTTl) together ensure that if the number and momentum sum rules 
(including the photon contribution) are satished for the PDFs in the proton, they are au¬ 
tomatically satished for the PDFs in the neutron^, regardless of the choices for and d^. 
Again, for our analysis, we choose uo and do to equal the initial up and down CT14 NLO va¬ 
lence distributions in the proton. Thus, from Eqs. 0-0 and flTTl) we can obtain the quark. 


^ This simple approximate isospin symmetry is broken by the inclusion of the elastic component of the 
photon PDF, since there is no corresponding elastic photon in the neutron. 
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gluon, and photon PDFs in both the proton and neutron in terms of the parametrization of 
the quark and gluon PDFs in the proton, plus the two additional parameters Au and Ad- 
For this paper, we shall make the further simplihcation of A^ = Ad, which corresponds 
to cutting off the radiation from both valence quarks at the same scale. With this choice, 
everything is specihed by one additional parameter, which can be taken to be the cutoff 
scale Qcut = Qu = Qd, dehned by A^ = Ad = In (QQ/Q^ut)- Alternatively, we can trade this 
parameter for the initial inelastic photon momentum fraction in the proton: 


Po = 


dxxff/piXyQo) 


( 12 ) 


For the remainder of this paper, unless otherwise specihed, the photon PDFs will be in this 
one-parameter radiation ansatz labeled by Pq. For comparison purposes, we will make one 
exception to this by dehning a “Current Mass” (CM) photon PDF, analogous to the MRST 
current mass PDF, and given by Ai = In (Qo/Qf) with = ruu = 6 MeV and Qd = md = 10 
MeV. For this choice the initial inelastic photon momentum fraction is determined to be 
Pq = 0.26%. For all other partons in our analysis we use the CT14 NLO initial distributions, 
except that to maintain a total momentum fraction of 1, we re-normalize the initial up, 
down, and strange sea-quark distributions, to account for the additional photon momentum 
fraction. Given that the relevant photon momentum fractions are very small, we hnd that 
this reduces the sea-quark distributions by typically less than 1%, and it is inconsequential 
in our analysis. (The sea-quark distributions are reduced by 0.9% for pg = 0.14%, and they 
are reduced by 1.6% for the CM photon PDF.) In Fig. [1] we plot the quantity xf^x^pp) 
for three representative photon PDFs, relative to the quark and gluon PDFs, at the scales 
Pf = 3.2 GeV and pp = 85 GeV. We note that the effect of the initial photon PDF and 
the QED evolution on the quark and gluon PDFs is imperceptible in these plots. For the 
photon PDF, we plot for pg = 0% (Qcut = Qo = 1.295 GeV) and for pg = 0.14% (Qcut = 71 
MeV), and for the CM photon PDF. 


III. CONSTRAINTS ON THE PHOTON PDF 

A. Constraints from the CT14 data set 

The constraints on the photon PDF from the DIS and Tevatron data, used in the CT14 
analysis, are relatively weak. These come from two main sources: isospin violation effects 
in nuclear scattering and constraints from the momentum sum rule. In general, isospin 
violation will arise through QED evolution, as well as from the initial conditions given by 
Eq. ffTTll . This isospin violation can be seen in Fig. O where we plot fd/n{x, pp) / fu/p{x, pp) 
and fu/n{x, pf)/ fd/p{x, pp) for several values of pp for the case where the initial photon PDF 
is zero, and for the case where the initial photon is the CM choice. Note that the isospin 
violation is small and most important at large x. Given that cuts of = Q^(l/x — 1) > 12 
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FIG. 1: Plots of xf{x,HF) for = 3.2 GeV (left) and = 85 GeV (right). Three representative 
photon PDFs are plotted: the “Current Mass” photon PDF (7 cm i red dotted), and photon PDFs 
with initial inelastic photon momenta fractions of pg = 0 and 0.14% ( 70 , blue dashed, and 70 . 14 , 
green dot-dashed, respectively). The effects of the different initial photon PDFs on the quark and 
gluon PDFs are imperceptible in these plots. 


GeV^, applied to enforce perturbativity in the calculations, typically require x < 0.2 — 0.4, 
we expect constraints from isospin violation to be small in the present data, as observed in 
the MRST analysis j^. 

Constraints from the momentum sum rule arise because any momentum carried by the 
photon implies less momentum available for the quark and gluon PDFs. In this way, con¬ 
straints from data on the colored parton PDFs indirectly impact the photon PDF. We have 


performed a preliminary analysis using the data sets included for CTIO 20|]. For a fixed ini¬ 


tial photon momentum fraction, with the photon PDF parametrized as discussed in Sec. II, 
we minimized the global by varying the quark and gluon PDFs. Using the usual CTEQ- 
TEA choice of < 100 tolerance, we obtain a limit on the photon momentum fraction 
of pg < 5.6% at the 90% confidence level, which is similar in magnitude to the results found 
by the MRST and NNPDF analyses. The best fit for the initial photon momentum fraction 
from this global analysis is pg = 1.2%, but with only a small change of = —7, relative 
to the fit with pg = 0%. For comparison, we find the elastic contribution to the initial 
photon momentum fraction, as calculated in the equivalent photon approximation, to be 

Unfortunately, this limit on the initial photon momentum fraction is much larger than one 
would expect for a photon PDF. In the analysis of the NNPDF group, additional constraints 
were made on the initial photon PDF by including LHC data on high-energy W, Z, and Drell- 














FIG. 2: Plots oi fd/n{x, fip) / fu/p{x, fJ-p) fJ-p)/fd/pix, fip) for fip = 1.3 GeV (black dots), 

3.2 GeV (red dashes), and 85 GeV (blue solid). The left plot is for zero initial photon momentum 
and the right plot is for the CM photon PDF. 

Yan production, and comparing with theoretical predictions that included photon initial- 
state contributions. Although the photon-induced contribution to these processes is small 
compared to the dominant quark-antiquark annihilation subprocess, the precision of these 
measurements was enough to substantially increase the constraints on the photon PDF j^. 
However, the small relative contribution of the photon-photon subprocess puts a stringent 
requirement on the precision needed for both experimental and theoretical analyses. Any 
small misjudgment of systematic errors on the experimental side, or uncalculated higher- 
order corrections on the theoretical side could have a significant effect on the extraction of 
the photon PDF. In particular, given that the initial photon PDF is nominally of order a, 
one might expect that the uncalculated 0{a‘^) quark-initiated contributions to Drell-Yan 
production would contribute at the same level as the photon-initiated contributions. For 
this reason, we consider a different experimental process, isolated photon production in DIS, 
to constrain the photon PDF. 


B. Calculation of the process ep ^ ej + X 

At the partonic level, the process of DIS with isolated photon occurs at LO through 
Compton scattering of a photon coming from the proton off the lepton, as shown in Fig.[3](a). 
Thus, this process probes the photon PDF at LO, having no large backgrounds with which 
to compete. However, the quark-initiated subprocess shown in Figs. [3](b) and (c), while 
formally suppressed by 0{a), is just as large because of the small size of relative to fq. 
In fact, if we consider the photon PDF to be 0(a), then the photon-initiated subprocess 
and the quark-initiated subprocess are actually the same order in a. Thus, the correct way 
to calculate the cross section for DIS with isolated photons is to include both subprocesses 
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(a) 




FIG. 3: Amplitudes for the process ep ^ 67 + X. For each diagram shown there is an additional 
diagram where the photon is emitted off the initial-state lepton or quark. 


consistently without double-counting. 

In the literature there have been two approaches to calculations of the process ep 
ey -f X. The calculation of MRST j^, which was preceded by studies of Bliimlein et 
al. |2ll-l23j. included just the photon-initiated contribution of Fig. |3]^a). The calculation 


of Gehrniann-De Ridder, Gehrmann, and Poulsen (GGP) 2J, [2^ included just the quark- 


initiated contributions of Figs. |3]^b,c). In the GGP analysis it was found convenient to make 
the Lorentz-invariant separation of the cross section into three components, depending on 
the fermion line off which the hnal-state photon is emitted: LL for emission off the lepton 
line, given by the square of the partonic amplitude in Fig. [3]^b); QQ for emission off the quark 
line, given by the square of the partonic amplitude in Fig. |3]|^c); and QL for the interference 
between the two sets of diagrams.'^ In the GGP calculation a cut on the outgoing quark 
was necessary to remove the divergence in the amplitude as the photon off-shellness went to 
zero in the LL amplitude. A hybrid calculation was also considered by the ZEUS Gollabo- 
ration in their analysis of the DIS-plus-isolated photon data 10|, where the LL component 
of the quark-initiated subprocess of GGP was replaced by the photon-initiated subprocess 
of MRST. 

In this section we introduce a consistent and systematic method of combining the photon- 
and quark-initiated subprocesses, which also reduces the factorization scale dependence of 
either calculation. First, consider the calculation of the differential cross section as a power 
series in a without consideration of the relative sizes of f.y and fq. It can be written as a 
convolution over partonic cross sections 


da = V / fa{^,PF)d&a , 
„ do 


(13) 


^ Note that each of the diagrams in Fig. factually represents two Feynman diagrams, where the final-state 
photon is emitted off the initial-state lepton or quark as well as off the final-state lepton or quark. 
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where each of the partonic cross sections can be expanded in a power series in a, 

da a = ^da^^^ . (14) 

n 

The diagram in Fig. [3](a) corresponds to a LO contribution (a = 7 ; n = 0), while the 
diagrams in Figs. [3](b,c) correspond to NLO contributions (a = g, g; n = 1). Through NLO 
in a the cross section can still be written as a sum of LL, QQ, and QL components, 


da = da^^^^ + da^^^^ + da^^^^ , 


where the LL component also includes the photon-initiated contribution, 



f'i{^,LF)da^ + fa{^,fiF)dai^^^ 


( 16 ) 


( 16 ) 


Using the modihed minimal subtraction (MS) scheme, we can factorize the initial-state 
singularity from the NLO quark-initiated subprocess into the dehnition of the photon PDF, 
leaving a NLO subprocess cross section. 


a 


aa^ aag -f- 


47r/i 




eF(l 


dzP^q{z) dd^^\zC) , 


(17) 


where the hrst term on the right-hand side, da^'^^\ is the hard partonic quark-induced 
subprocess, and the second term is the collinear subtraction counterterm. Here, we have 
distinguished the renormalization scale d-R from the factorization scale d'F-, and we explicitly 
note that the initial-state collinear singularity cancels within the LL component. Calculating 
everything in dimensional regularization with 4—2e dimensions, the combined LL component 
of the NLO quark-initiated subprocess cross section, is finite as e —)■ 0. 

In principle, there are additional virtual and real contributions at NLO in a, besides 
the quark-initiated contributions. However, all other NLO terms are proportional to f^, 
which is in fact suppressed by an amount of order a relative to fq, as seen in the previous 
section. Thus, by keeping the photon-initiated contribution at LO and the quark-initiated 
contributions at NLO, and including the collinear-subtraction counterterm of Eq. flTT)) we 
have a well-defined calculation in the MS scheme, while maintainin g all contributions of the 
same size in a. Note that this is reminiscent of the ACOT scheme 26(] for including charm 
and bottom quark PDFs, although for the photon, there is no equivalent of the fixed-flavor 
scheme, due to its zero mass and the consequently nonperturbative photon PDF. 

The calculations of the QQ and QL components are identical to those in the GGP calcu¬ 
lation. For the kinematics of interest to us, the QL component is negligible, but it is included 
for completeness. The QQ contribution has a final-state singularity, when the photon and 
final-state quark become collinear, which can be handled by including a fragmentation con¬ 
tribution and an associated collinear subtraction counterterm in the MS scheme. Thus, we 
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obtain 




27 r I J er(l — e) 


dzP'yqi^Z^ dCeq—^e'Q' 


Q'=k'lz 


dz ID'yjQi (^Z^ /^/) d&eq—^e'Q' 


Q'=k'lz 


( 18 ) 


where is the hard partonic subprocess, ddeq^e'Q' is the LO subprocess cross section 

for eq —)■ e'Q', and D^/qi{z, ^p) is the fragmentation function for finding the photon (with 
momentum k') in the quark Q' (with momentum Q'), with momentum fraction z at the 
fragmentation scale qif. The singularities in the hrst two terms on the right-hand side of this 
equation cancel as e —)■ 0. We will discuss the choice for the fragmentation function in the 
next subsection. 

We have calculated the LL and the QQ components of the cross section, using the sub¬ 
traction method to handle the collinear divergences 27| in the hard cross-section term. 
Following this method, we subtract a term with a two-particle final state mapped onto a 
3-particle final-state, with the third particle phase space unintegrated, and then add back the 
exact same term with the third particle phase space integrated out. The subtraction term 
is designed to have the same collinear-singular limit as the hard term in the same region of 
phase space of the third particle, so that the hard cross section term minus the subtraction 
is integrable in d = 4 dimensions, while the 1/e singularities in the remaining terms cancel. 
Using this method we obtain for the LL quark-initiated contribution: 






da 


il,LL) 

q'(sub) 


4- dd 


(1,LL) 

g(AP) 


(19) 


where 


“‘^g(AP) 


a 


dz — 
27r 


P^q{z)\n 


g^(l - zf 


+ 6 ^^ ddf(zO- 


( 20 ) 


The hard term minus the subtraction term can be written as an integral over the phase space 
of the additional quark: 


aaq “<^g(sub) I “ 




f27r 


dcj) 


if ( dcr'''**’ 


IP"I S''"- 


( 21 ) 




dw dcj) 


where ta = (1 — cosd)/2 with 6 and 0 the scattering angles of the final-state quark in the 
initial parton-parton center-of-momentum frame. Note that the hard term on the right-hand 
side is treated with three-body final-state phase space, while the subtraction term is treated 
with two-body final-state phase space. 
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Similarly, we obtain for the QQ contribution: 




(1,QQ) 

9 (frag) ’ 


where 


( 22 ) 


da 


(1,QQ) 

9 (frag) 


dzl D^/Q,(z,izf) 


(23) 


a 


P^g{z) In ( I + 


The hard term minus the subtraction term can be written 


da. 


eq^e'Q' 


“%(sub) ' “ 


dw I — 271——- 

2tt \ dw d(j) 


a 


dz 2yj-.j2) ddeq^e'Q' 


(24) 


where in this case we have found it convenient to use a different parametrization of the 
hnal-state quark phase space. Letting qe —)■ g'e'y be the hard partonic subprocess, we use 
tc = (1 — cos6')/2, where 9 is the angle between q' and 7 in the q'e' center-of-momentum 
frame and (j) is the azimuthal angle between the qej plane and the q'Pj plane. As before, 
the hard term on the right-hand side is treated with three-body hnal-state phase space, 
while in this case the subtraction term and the fragmentation term are treated in the limit 
where the hnal-state photon and quark are collinear, with momenta satisfying k' = zQ' and 
g' = (1 - z)Q'. 


C. ZEUS Experimental Cuts and Photon Isolation 


The ZEUS experiment 10|] used proton and lepton beam energies of Up = 920 GeV and 
El = 27.5 GeV, respectively, corresponding to a center-of-mass energy and rapidity of 


\fs = 2 EpEi = 318 GeV 

K = ilti^ = 1.76, (25) 

respectively. (We neglect the proton mass, mp, in all calculations here.) For the process 
ep —)■ ey -|- X, with momentum satisfying, ^ + p = P + k' + p'x, one can dehne the standard 
DIS variables that describe the kinematics of the scattered lepton, = — (£ — U)^, y = 
p • (£ — U)/(p • £), and X = / (sp). The ZEUS Gollaboration measured distributions for two 

leptonic variables Q and x, and for two photonic variables, Ej^^ and the transverse energy 
and pseudorapidity of the photon, respectively. The collaboration combined data that were 
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59.1% e~p and 40.9% e+p scattering. Note that the sign of the charged lepton has no effect 
on the LL or QQ components of the calculation, but the combination of the two charged 
lepton contributions produces a signihcant cancellation of the already small QL component, 
so that it is negligible in the analysis. 

The kinematic region dehned by the experiment was 


10 GeV^ 

< 


< 

350 GeV 



A 

> 

10 GeV 

139.8° 

< 

O'e 

< 

171.8° 

4 GeV 

< 


< 

15 GeV 

1 

o 

< 

Vi 

< 

0.9 . 


(26) 


The cut on the hnal-state lepton angle, 0^, can be written in terms of its rapidity as 


- 2.6355 < < -1.0053 . 


(27) 


There are two additional cuts that require discussion. The experimentalists reported a cut 
of Wx > 5 GeV, where Wx = + l — — Naively, this cut looks problematic because it 

would remove the photon-initiated contribution, which occurs at exactly Wx = 0. However, 
upon closer investigation it appears that this cut was only applied to the theoretical and 


Monte Carlo calculations. To quote from Ref. 28| , “ The keen reader will note that no such 


cut was applied at the detector level. This proved impossible due to the poor Wx resolution 
at detector level and poor description of the data by MC...”. The relevant detector-level 
cut was the requirement of at least one reconstructed track, well separated from the lepton, 
which was used to ensure some hadronic activity and to remove deeply virtual Compton 
scattering events. After this cut, it was found that the number of events in the Monte Carlo 
calculations with Wx < 5 GeV was negligible. For our purposes, we interpret the forward 
track cut as a requirement to tag inelastic events, and we include no explicit Wx cut. 

Note that the forward track cut should equally remove elastic isolated photon events, 
and so remove the contribution from the elastic part of the photon PDF. In this way, the 
ZEUS data probe only the inelastic part of the photon PDF, and therefore, we only include 
this inelastic contribution in comparison with the experimental data. In doing so we have 
made the approximation /-y, inclusive (a;, Q) ~ A, elastic (ai, Q) + A, inelastic (ir, Q); he., the elastic 
and inelastic components of the photon PDF evolve separately. This approximation is good 
because A,eiastic(a^) Q) changes very little from Qq to Q due to the rapid falloff of the proton 
electromagnetic form factor, while the inelastic contribution evolves additively. 


Ai inelastic (ai, Q) ~ A , inelastic ( 


ai,Qo)+ y, / 

i JQI 


dQ^ a 


p2 p(0) 


° (A + A-i)(uQ), (28) 


up to corrections suppressed by extra factors of a. We have verihed by explicit calculation 
that this additive approximation replicates the consistently evolved inclusive photon PDF, 
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with errors that are far smaller than other theoretical uncertainties that we will discuss 
below. 

The second additional important cut is the isolation cut on the photon, enforcing that 
90% of the energy in the jet containing the photon belongs to the photon, where jets are 
formed with the cluster algorithm with parameter R = 1.0. We will model this isolation 
cut in our calculation in two different ways. First, we can model the experimental cut 
at the parton level, requiring E^/{Eg + E^) > 0.9 if the photon-quark separation satishes 
r = = (^{Arj^qY + {A(j)^q)^^ < 1. For later reference, we call this the “sharp” 

isolation cut. Since this does not completely remove the quark-photon collinear singularity, 
the theoretical calculation of the QQ component will depend on the choice of the quark-to- 


photon fragmentation function D^/qi{z, fif 
determined by the ALEPH Collaboration 


For this, we use the LO fragmentation function 


29|, parametrized by 




a 


P^q{z) 


In 




where /Uq = 0.14 GeV and Cq = —13.26. Note that this parametrization of D^/qi{z, ^f) 
is an exact solution to the evolution equation at 0(a), so that the dependence on the 


fragmentation scale /i/ cancels exactly in our calculation. In the GGP analysis 2^ , other 


parametrizations of D-y/Qi{z, fif) with different assumptions were compared, with only a small 
effect on the calculated cross sections. 

One of the disadvantages of having a dependence on the quark-to-photon fragmentation 
function, in addition to the uncertainties due to the phenomenological £t to the function, 
is that it assumes that the cross section is inclusive in the fragmentation remnant. In 
our calculation, the combination of the experimental constraints on the photon and on the 
lepton indirectly imposes constraints on the remnant quark in the process. Therefore, we 
also consider an alternative model of the experimental isolation cut, replacing it with a 
“smooth” isolation cut 30|, so as to avoid the necessity of the fragmentation contribution. 
The smooth isolation cut is given by requiring that the hadronic energy Eh inside all cones 
of radius r < R around the photon direction satisfy 


Eh < eE^ 


1 — cos r 
1 — cos R 


(29) 


where we take R = 1 and e = 1/9. These values of R and e are chosen to ensure that 
the photon will contain at least 90% of the energy inside a cone of r = 1.0 centered on 
the photon, just as for the experimental isolation cut. However, the smooth cut does not 
translate exactly to the experimental isolation cut, because it requires the photon to carry a 
greater fraction of the energy as the cone size r becomes smaller. In practice, because of this, 
the theoretical calculation with the smooth cut is better behaved than with the sharp cut. In 
addition, the smooth isolation cut is more restrictive than the sharp isolation; for a strictly 
positive-definite differential cross section, the smooth isolation prescription must always give 
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FIG. 4: Differential distributions for a zero initial inelastic photon PDF, using the smooth isolation 
prescription. The various bands display a variation in factorization scale between 0.5F^_L..y < ^ 

2E_i_'y and correspond to the total prediction (light gray solid), the QQ component (blue hashed), 
the LL component (dark red solid), and the photon-initiated contribution only (green hashed). 
Also shown are the ZEUS data points with combined statistical and systematic errors. 

a smaller predicted cross section. In this way, a comparison of the two calculations can give 
some indication of the theoretical uncertainty due to the isolation cut. 

Finally, we note that the jet-clustering algorithm also includes the electron, so that the 
isolation cut effectively imposes AR^e > 1. 


D. Comparison with the data and constraints on the photon PDF 

Before discussing the comparison of the theory with the data, it is useful to understand the 
theoretical uncertainties of the calculation by studying the factorization scale dependence 
and the dependence on the isolation prescription. In Fig. 0] we plot the differential cross 
sections for da/dE±^ and da/dr]^ as functions of Ej_'y and respectively, while varying 
the factorization scale within 0.5F'_L.y < fip A 2E±^. Here we have used the sharp isolation 
cut and calculated with zero initial inelastic photon PDF at Qo = 1.295 GeV. The four 
bands on each of the two plots correspond to the photon-initiated contribution only (green 
hashed band), LL component (dark red solid band), QQ component (blue hashed band), 
and the total calculation (light gray solid band). The QL contribution is imperceptible on 
the scale of these plots. From these plots we learn several important facts. First, the scale 
dependence of the LL component is reduced dramatically compared to the photon-initiated 
contribution alone. This large scale dependence of the photon-initiated contribution cancels 
greatly with that of the collinear subtraction counterterm in the combined LL component. 
Second, the LL and QQ components dominate in different regions of phase space. For 
instance, the cross sections are most sensitive to the LL component, and consequently to 
the photon PDF, at large E±^ and small rj^. Thus, the shapes of these distributions can give 
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FIG. 5: Differential distributions for a zero initial inelastic photon PDF with the factorization 
scale varied between 0.5Fi_L.y < < 2E±^. The blue hashed band is calculated using the smooth 

isolation prescription and the red hashed band is calculated using the sharp isolation prescription. 
Also shown are the ZEUS data points with combined statistical and systematic errors. 

information about the nonperturbative contribution to the photon PDF. Finally, we note 
that the scale dependence of the QQ component is still large, being only LO in a*, and it 
dominates the overall scale uncertainty of the theoretical calculation. 

In Fig. [5] we plot the total predictions of the same two distributions, again for zero ini¬ 
tial inelastic photon PDF, but now comparing the two different isolation prescriptions in the 
theoretical calculation. In these plots we show the predictions with the smooth isolation pre¬ 
scription (blue hashed band) and the sharp isolation prescription (red hashed band), again 
varying the factorization scale within 0.5E±^ < < 2E±^. The hrst thing to note here is 

that the difference between the predictions is about the same size as the scale uncertainty, 
with similar dependence on the kinematic variables. Another striking feature is that the pre¬ 
diction using the smooth isolation prescription is uniformly larger than that using the sharp 
isolation, in contrast to expectations. This is probably due to incomplete cancellations in 
the sharp isolation calculation between the large negative collinear fragmentation contribu¬ 
tion and the positive real contribution, due to indirect constraints on the emitted hnal-state 
quark in the real emission contribution. Presumably higher-order QCD corrections will af¬ 
fect the predictions for both isolation predictions, to resolve this issue. As noted previously, 
our calculation is only LO in ag] we expect both the factorization scale uncertainty and the 
isolation prescription discrepancy to be reduced at NLO in In any event, we will use the 
two isolation prescriptions, as well as the factorization scale dependence, as a measure of the 
theoretical uncertainty of our calculation. 

With this understanding of the theoretical uncertainties of the calculation we can now 
compare the ZEUS data against predictions for the differential distributions, while varying 
the initial inelastic photon momentum fraction pg of the photon PDF, described in Sec. IIIBI 
For this analysis, the initial quark and gluon PDFs are just the CT14NLO PDFs, except that 
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FIG. 6: Differential distributions in the photon variables, and with the smooth isolation 
prescription, with factorization scale /tf = 0.5iil_L-y. The curves, from bottom to top are with initial 
inelastic photon momentum fractions of = 0% (black dashed), 0.1% (blue dotted), 0.2% (red 
dot-dashed), and for the CM photon (green solid). Also shown are the ZEUS data points with 
combined statistical and systematic errors. 

the sea-quark normalizations are rescaled in order to maintain a total momentum fraction 
of 1. This rescaling has little effect in our analysis, because the photon momentum fractions 
considered here are small. For instance, a photon momentum fraction of = 0.14% induces 
a reduction of the sea-quark momentum by only 0.9%, while the CM photon PDF induces 
a reduction of the sea-quark momentum by 1.6%. At this stage of the analysis we have not 
reht the quark and gluon PDFs, since the ep ^ ep -\- X process is dominantly sensitive to 
the photon PDF directly, whereas the indirect sensitivity through changes in the quark and 
gluon PDFs is negligible. 

In Fig. |6]we plot the differential cross sections for da/dE±^ and da/dp^ as a function of the 
photon variables, and using the smooth isolation prescription, with a factorization 
scale of pf = 0.5i?F7- The curves, from bottom to top are with initial inelastic photon 
momentum fractions of pg = 0% (black dashed), 0.1% (blue dotted), 0.2% (red dot-dashed), 
and for the CM initial photon (green solid), which has initial momentum fraction 0.26%. 
Also shown are the ZEUS data points with combined statistical and systematic errors. With 
these choices of /if and the isolation prescription, we see that the theory can £t the data 
well for Pq 0.1%. On the other hand, the theory hts poorly for the CM initial photon, 
overshooting the data at large Ex_^ and small p^. Of course, the best fit for pg is correlated 
with the choice oi pp and the isolation prescription. However, since these choices tend to 
move the curves up or down uniformly, it is still possible to constrain the initial photon 
PDF by the shape of the distributions. In particular, it is impossible to get a good £t to the 
prediction using the CM initial photon PDF, regardless of the choices oi pp and the isolation 
prescription. 

In Fig. [7] we plot the differential cross sections for da/dQ'^ and da/dx as a function of the 
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FIG. 7: Differential distributions in the lepton variables, and x, with the smooth isolation 
prescription, with factorization scale fip = O.SQ. The curves, from bottom to top are with initial 
inelastic photon momentum fractions of = 0% (black dashed), 0.1% (blue dotted), 0.2% (red 
dot-dashed), and for the CM photon (green solid). Also shown are the ZEUS data points with 
combined statistical and systematic errors. 

lepton variables, and x against the ZEUS data, using the exact same theoretical choices 
and initial inelastic photon PDFs as for the previous plot. In this case we see that it is 
impossible to fit the data, regardless of the initial photon PDF or the choices of factorization 
scale and isolation prescription. In particular the theory fits the data very poorly at small 
X and Q^. In fact, we note that prediction for the smallest bin in x is far from the ZEUS 
data point, and is essentially independent of the initial photon PDF. We expect that the 
predictions in these bins are highly sensitive to higher-order QCD radiation, so that it is 
difficult to £t the full lepton distributions with a hxed-order calculation. 

The fact that these fixed-order calculations are more reliable for the photon distributions 
than for the lepton distributions can be seen further by looking at the phase-space constraints 
for the two sets of variables. In Fig. Owe show plots of the constraints on the photon variables 
and and on the lepton variables and x. In these hgures, the dashed lines indicate 
the bins that are plotted by the ZEUS data. The combined dark red-l-light blue regions 
indicate the regions of phase space allowed by the ZEUS kinematic constraints of Eq. fl27|) . 
for the fully inclusive event, whereas the dark red region only is allowed for the LO photon- 
initiated subprocess. For the photon distributions, the constraints are dominated by the 
photon cuts on and 7]^, with only a small cut on the photon-initiated contribution in 
the upper-left corner due to the requirement of < —1.0053. Thus, all of the bins in E_i^ 
and 7]^ have a large photon-initiated contribution and can be considered very inclusive. In 
contrast, for the lepton distributions the additional photon constraints have a large effect in 
many of the bins. For instance, the photon-initiated contribution to the smallest bin is 
largely removed by the requirement of Ey^^ > 4 GeV, and the photon-initiated contribution 
to the smallest x bin is completely removed by the requirement of > —0.7. These bins are 
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FIG. 8: Kinematic limits on and r}^ (left) and and x (right). The dashed lines indicate the 
bins that are plotted by the ZEUS experiment. The dark red region only is kinematically accessi¬ 
ble by the photon-initiated contribution, while the dark red-|-light blue regions are kinematically 
allowed in general. 

dominated by events with additional particles in order to satisfy the kinematics, so we would 
not expect our hxed-order calculation to do well at predicting the lepton distributions. 

Based on these arguments we will use only the distributions in the photon variables E_i_^ 
and to constrain the initial inelastic photon PDF, for a total of eight data points. We also 
reiterate that the constraints due to the remaining CT14 experimental data set are much 
weaker than these data, and are neglected in the present analysis. We dehne the chi-squared 
function for these data points by 



where T,, Di, and cTj are the theory prediction, the experimental measurement, and its 
combined statistical and systematic error for the data point i. In Fig. [9] we plot versus 
the initial inelastic photon momentum fraction pg for both the smooth and sharp isolation 
prescriptions and for several values of the factorization scale fip. Note that the value of pg 
determined by the minimum of depends signihcantly on the isolation prescription and on 
the factorization scale, giving best hts for the initial momentum fraction that can vary from 
less than 0 to above 0.1%. In addition, due to this theoretical uncertainty in the current 
calculation, it is not possible to unambiguously determine an error band on pg, using the 
standard CT approach of applying some tolerance criterion on the rise in the x^ around 
the best £t. However, from Fig. [9] we do see that not all choices of theoretical parameters 
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FIG. 9: Plots of ^ versus initial inelastic photon momentum fraction using the smooth isolation 
prescription (left) and the sharp isolation prescription (right) for factorization scales pp = 

0.5E±ry, and 0.35ill_L.y. The horizontal line at = 13.36 is the 90% confidence level limit for 
eight data points. 

are able to fit the shape of the data points equally well. Therefore, we can determine a 
conservative limit on the value of Pq by requiring that the data and theory not disagree 
beyond some level. A distribution with eight data points will have < 13.36 at the 90% 
conhdence level.^ Therefore, we define that any theoretical prediction with > 13.36 is 
ruled out as a bad £t to the data at the 90% confidence level. It is impossible to satisfy this 
criterion for > 0.14% for either choice of the isolation prescription and for any value of 
pp. Furthermore we find that the CM choice of the photon PDF has > 46 for any choice 
of isolation and factorization scale and so is ruled out by this data. 

Thus, we find our maximal initial inelastic photon PDF to have Pq = 0.14% at the 90% 
conhdence level. Of course, the exact value of the momentum fraction is correlated with the 
shape of the initial photon PDF. From Fig. |H]we see that the ZEUS DIS-plus-isolated-photon 
data constrains the photon PDF in the kinematic region given roughly by 10“^ < x < 2- 10“^ 
for 16 < < 300 GeV^. Outside of this region the photon PDF is very weakly constrained, 

but we believe that the radiative ansatz gives a reasonable expectation for its overall shape. 
As for the minimal possible value of the initial inelastic photon momentum fraction, it could, 
in principle, be negative, which is not ruled out by the analysis of this section. For instance, 
one could begin the evolution with zero initial photon PDF at a lower value of the scale Qo- 
However, we take the reasonable assumption that it should be nonzero at the low scale of 
Qo = 1-295 GeV. We thus propose the initial PDFs with p^ = 0% and p^ = 0.14% as our 
90% G.L. photon PDFs. A similar analysis gives p^ < 0.11% at the 68% conhdence level, 
but the data are still consistent with p^ = 0% at the 68% G.L. 


® As a comparison, the change in the total for the remaining 2947 data points used in the CT14 analysis 
is Ay^ = —2.3 in going from pg = 0% to pg = 0.14%. 
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FIG. 10: Comparison of various NLO photon PDFs at the scale Q = 3.2 GeV: CT14QED with 
Pq = 0% (solid green), CT14QED with pg = 0.14% (short-dashed black), MRST2004QED0 us¬ 
ing current quark masses (dotted orange), MRST2004QED1 using constituent quark masses (dot- 
dashed brown), and NNPDF2.3QED with ag = 0.118 and average photon (long-dashed blue). 


In Fig. do] we compare, at the scale Q = 3.2 GeV, the CT14QED photon PDFs with 
Po = 0% and pg = 0.14% against the NLO MRST2004QED photon PDFs, using the current 
quark masses (labeled MRSTO) and using the constituent quark masses (labeled MRSTl), 
and against the NLO NNPDF2.3QED average photon PDF with = 0.118. We should 
emphasize that the CT14QED photon PDFs only contain the inelastic contribution in these 
plots. The NNPDF2.3 average photon PDF has more structure in its shape at large and 
small values of x than do the other PDFs, but it is still consistent with the ZEUS data in 
the X range that is probed by the experiment. 

In Fig. dllwe compare the same set of photon PDFs at the higher scale of 85 GeV and the 
very high scale of 1 TeV. Here we can make some very interesting observations. The most 
obvious thing in these hgures is that the GT and MRST photon PDFs become very similar at 
large whereas the NNPDF photon PDF is distinctly different and much smaller at small 
values of x. This difference is due to the different approaches to the evolution of the PDFs 
taken by the different groups. Whereas in the MRST and GTEQ-TEA approaches, the QGD 
and QED scales are chosen to be identical and evolved together, in the NNPDF2.3QED PDFs 


that are included in the LHAPDF library 3l|, the QGD and QED scales are separate and 


the two scales are evolved successively; hrst the QED scale is evolved from Qq to Q and then 
the QGD scale is evolved from Qq to Q. As discussed in Ref. jsl, the successive evolution of 
QED and QGD differs from the combined evolution by terms that are subleading by 0{aas) 
and can induce large unresummed logarithms between the two scales. This difference in the 
evolution at small x is also seen to be consistent with the behavior seen in the right panel of 
Fig. 2 in Ref. j^, where the NNPDF photon PDF also is smaller at small x and large than 


when it is evolved using the code partonevolution [7], Il8|. We expect that the difference 
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FIG. 11: Comparison of various NLO photon PDFs at the scales Q = 85 GeV (left) and Q = 1 TeV 
(right): CT14QED with pg = 0% (solid green), CT14QED with pg = 0.14% (short-dashed black), 
MRST2004QED0 using current quark masses (dotted orange), MRST2004QED1 using constituent 
quark masses (dot-dashed brown), and NNPDE2.3QED with Og = 0.118 and average photon (long- 
dashed blue). 

between the NNPDF2.3QED photon PDF and the other photon PDFs at high Q would be 
less significant if the NNPDF2.3QED PDFs were evolved from the low scale simultaneously 
in QED and QCD. 

Another observation from Fig. [HI concerning the CT14QED and MRST2004QED photon 
PDFs is that the impact of the initial photon distribution becomes less significant as 
increases and more photons are produced through radiation off the quarks. From these 
plots we see that the fractional deviation between the different photon PDFs decreases 
with increasing Q^. In fact at very small x and large the differences in the sea-quark 
distributions of the PDFs presumably have more impact on the photon PDF than does the 
initial photon distribution. 


IV. DISCUSSION AND CONCLUSION 

In this paper, we have presented CT14QED, which is the hrst set of CT14 parton distri¬ 
bution functions obtained by including QED evolution at leading order (LO) with next-to- 
leading-order (NLO) QCD evolution in the global analysis by the CTEQ-TEA group. This 
development will provide better theory predictions to compare with the precision data, such 
as Drell-Yan pair production, measured at the LHC. The CT14QED PDFs are based on the 
CT14 NLO initial distributions with the addition of an initial photon PDF. (There is also 
an inconsequential rescaling of the quark sea PDFs, in order to maintain the momentum 
sum rule.) The inelastic contribution to the photon PDF is parametrized at the initial scale 
Qo using a generalization of the radiative ansatz introduced by the MRST group in their 
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previous study. The initial photon PDF then depends on two independent parameters [cf. 
Eq. ([I])], which are related to the scales at which the radiation off the up and down valence 
quarks is cut off. However, given the weak constraints on the initial photon PDF we hnd it 
convenient at this time to set the scales equal, so that the initial photon PDF is parametrized 
by a single parameter, which we take to be the momentum fraction carried by the inelastic 
photon at the initial scale Qq. For comparison purposes we have also dehned a “Current 
Mass” (CM) photon PDF, comparable to the MRST current mass PDF, for which the initial 
photon momentum fraction is = 0.26%. 

A set of neutron PDFs can also be obtained with a small amount of isospin breaking, 
suggested by the radiative ansatz applied to hrst order in a, and which automatically ensures 
that the number and momentum sum rules are satished. However, as previously seen by 
both the MRST and NNPDF groups, we hnd that the constraints from isospin violation 
effects (generally small and most important at large x) in nuclear scattering and from the 
requirement of the momentum sum rule, imposed by the DIS and Tevatron data in the CT14, 
are relatively weak. 

Thus, in order to constrain the photon PDF, we focused on the scattering process ep —?■ 
ejX, which was measured by the ZEUS experiment at HERA. This process is dominantly 
sensitive to the inelastic photon PDF directly, with negligible indirect sensitivity through the 
modihcation of the quark and gluon PDFs by QED effects. It also has the advantage that 
the initial-state photon subprocess contribution occurs at leading order, so that it does not 
compete with other much larger contributions. In this paper we have produced for the hrst 
time a consistent and systematic calculation for this cross section that combines both the 
photon- and quark-initiated subprocesses, and simultaneously reduces the factorization scale 
dependence of either calculation. Details of this calculation were presented in Sec. IHIB1 The 
photon isolation cut, which required that the hnal-state photon must contain at least 90% of 
the energy in the jet to which it belongs (where jets are formed with the fcr cluster algorithm 
with parameter R = 1.0), was modeled using two different models of photon isolation. We 
used the two different isolation prescriptions, as well as the factorization scale dependence 
as a measure of the theoretical uncertainty of our calculation. 

By comparing the ZEUS data for the distributions of transverse energy and pseudorapidity 
of the hnal-state photon against our calculation of the differential distributions, we were able 
to constrain the initial inelastic photon momentum fraction inside the proton to be < 
0.14% at the 90% conhdence level. Hence, the CM choice of photon PDF has been ruled out 
by this data. For completeness, we also compared the CT14QED PDFs to some of the NLO 
(in a^) photon PDFs published by the MRST and the NNPDF groups. Phenomenological 
applications of the CT14QED PDFs will be discussed in future publications. 

As shown in Fig. |5l the theoretical uncertainties due to the factorization scale dependence 
and the isolation prescription are currently larger than the experimental uncertainties of the 
Zeus data. Thus, extending our calculation to NLO in should be able to further constrain 
the initial photon PDF. This is a project that we are currently undertaking. 
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Parametrizations for the (inelastic photon) CT14QED PDF sets (both proton and neutron 


versions) will be distributed in a standalone form via the CTEQ-TEA Web site 32 


or as 


a part of the LHAPDF6 library [Mf- For backward compatibility with version 5.9.X of 
LHAPDF, our Web site also provides CT14 grids in the LHAPDF5 format, as well as an 
update for the CTEQ-TEA module of the LHAPDF5 library, which must be included during 
compilation to support calls of all eigenvector sets included with CT14 |33|. We will also 


distribute sets with the inclusive photon PDFs, CT14QEDinc. For the proton, CT14QEDinc 
at the initial scale Qo is the sum of the (inelastic) CT14QED and the elastic component of 
the photon PDF, given by the Equivalent Photon Approximation. The proton CT14QEDinc 
PDFs are then evolved from Qq to Q as discussed in Sec. Ill Al For the neutron, CT14QEDinc 
is equal to CT14QED, since the neutron has zero electric charge, and therefore it has no 
elastic component of the photon PDF. 
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Appendix: Comparison of our QCD plus QED evolution code with other codes 


We have checked our code against the public evolution code partonevolution 00, 
which also solves the evolution equations at LO in QED and NLO in QCD. The main 
difference (other than technicalities of implementation) in the partonevolution code is 
that it also includes the charged leptons as partons in the proton. Using the toy model 
of Ref. 3^ with zero initial photon PDF, ensuring all input parameters agree, using the 
same formulation for the running of ag and a, and removing the lepton PDFs and their 
contribution to from partonevolution, we find excellent agreement between the two 
programs. Evolving from Qq = 1-3 GeV to fip = 100 GeV, we find differences of less than 
0.2% over most of the range of x for all of the PDFs, including the photon. Reinstating the 
lepton PDF contribution to the evolution equations in partonevolution, which in principle 
should be included for consistency, we find their effects on the quark and gluon PDFs to 
be negligible, with changes of less than 10“^%. The only noticeable effect is the reduction 
of the photon PDF by about 1% with mild x dependence, due to splitting of the photons 
into lepton-antilepton pairs. This is presumably much less than other uncertainties in our 
analysis, so it is reasonable to leave out the lepton PDF contribution in our code. The 
percent difference in our prediction for the photon PDF, relative to the partonevolution 
prediction with and without the inclusion of the lepton PDFs, is shown in Fig. [T2j 

We have also checked our code against the program APFEL jsj, which includes QED at LO 
and QGD at up to NNLO. The main difference in the APFEL program is that the QGD and 
QED factorization scales can be taken to be independent, and the evolution with respect to 
each scale can be done successively. However, using the setting “QavDP” in APFEL averages the 
two possible orderings for performing the evolutions, which should agree with our approach 
to O(a^). In addition, the APFEL code has been recently updated to allow the simultaneous 
evolution of the QED and QGD scales, using the “QUniD” setting. We have compared our 
code with APFEL 2.6.0, starting with the GTIONLO PDFs 20| with zero initial photon PDF 
at Qq = 1.3 GeV and evolving in QED at LO and QGD at NLO to fip = 100 GeV. We 
have done the comparison using both the “QavDP” and the “QUniD” settings for APFEL. We 
obtain excellent agreement for the quark and gluon PDFs, with differences of less than 0.2% 
over most of the range of x for both APFEL settings. In Fig. [13] we show the results for the 
photon PDF. We obtain pretty good agreement with APFEL with the “QavDP” setting, with 
differences of less than 2.5% except at large x > 0.1. This is consistent with the 
differences expected in the different evolution procedures. We obtain excellent agreement 
with APFEL with the “QUniD” setting, with differences of less than 0.34% over the full range 
of a; > 10“® shown. We note that in Fig. [13] we replace the evolution subroutine for a in 
the APFEL program with the code used in the GT global analysis code; however, using the 
original a subroutine in APFEL still gives differences of less than 1% for the evolved photon 
PDF over the full range of x when using the “QUniD” setting. This is certainly much smaller 
than the uncertainties in the initial photon PDF itself. 
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FIG. 12: Percent difference in our prediction for the photon PDF, relative to the partonevolution 
prediction, for the toy model parametrizations evolved from Qq = 1.3 GeV to fip = 100 GeV, 
as discussed in the text. The solid red curve uses the partonevolution calculation with lep¬ 
ton PDFs removed from the evolution, while the dashed blue curve includes the lepton PDFs in 
partonevolution. 
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FIG. 13: Percent difference in our prediction for the photon PDF, relative to the prediction by 
APFEL with the “QUniD” setting (red solid) and “QavDP” setting (blue dashed). The PDFs are 
CTIONLO with zero initial photon PDF, evolved from Qq = 1.3 GeV to fip = 100 GeV. 
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We have not checked our code directly against the MRST evolution code or the recently 
developed QCDNUM+QED evolution code but we do note that the comparison between 


QCDNUM+QED and APFEL “QavDP” for the evolution of the photon PDF in Ref. 


35 


looks 


qualitatively similar to the results that we have found in Fi g. 11^ In addition benchmarking 
studies between APFEL and these two evolution codes in Ref. [3^ show agreement at a similar 
level to that which we have found with our code here. 
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